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Abstract 

We study the flavour singlet pseudoscalar mesons from first principles 
£^ , using lattice QCD. We explore the quark content of the r\ and r\ mesons 

■ and we discuss their decay constants. 

o : 

(N ■ 

1 1 Introduction 

\o : 

There is considerable interest in understanding hadronic decays involving 77 and 
77' in the final state. The phenomenological study of hadronic processes involving 
flavour singlet pseudoscalar mesons makes assumptions about their composition. 
^ ' Here we address the issue of the nature of the 77 and 77' from QCD directly, 

'-H • making use of lattice techniques. 

O 1' Lattice QCD directly provides a bridge between the underlying quark de- 

scription and the non-perturbative hadrons observed in experiment. The am- 
plitudes to create a given meson from the vacuum with a particular operator 
made from quark fields are measurable, an example being the determination of 
f„. It also allows a quantitative study of the disconnected quark contributions 
that arise in the flavour singlet sector. The lattice approach provides other in- 
formation such as that obtained by varying the number of quark flavours and 
their masses. 

In the case of pseudoscalar mesons, the chiral perturbation theory approach 
also provides links between a quark description and the hadronic states. For the 
pion, this has the well known consequence that the decay constant /„- describes 
quantitatively both the /iv and 77 decays. For the flavour singlet states (77 and 
rf), the situation is more complicated The axial anomaly now involves a 
gluonic component and the definition of decay constants is not straightforward. 
From the chiral perturbation theory description, one expects the mixing of 77 
and 77' to be most simply described in a quark model basis. In the flavour singlet 
sector, for pseudoscalar mesons, we then have contributions to the mass squared 
matrix with quark model content (uu + dd) / \/2 and ss (which we label as nn 
and ss respectively): 
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V^2~2-ns ^^ss -Ess 

Here m corresponds to the mass of the flavour non-singlet eigenstate and 
is the contribution to the mass coming from connected fermion diagrams while 
x corresponds to the contribution from disconnected fermion diagrams. In the 
limit of no mixing (all x — 0, the OZI suppressed case), then we have the 
quenched QCD result that the rj is degenerate with the 7r meson and the 77' 
would correspond to the ss pseudoscalar meson. This is not the case, of course, 
and the mixing contributions x are important. 

Using as input m nni m SSl and m^i , the three mixing parameters a; cannot 
be fully determined. It is usual to express the resulting one parameter freedom 
in terms of a mixing angle, here defined by 

V = r) nn cos 4> - r/ ss sin <j> r\ = r\ nn sin 4> + r/ ss cos (f> (2) 

We show the resulting values of the mixing parameters x in the figure (the input 
value for m ss will be discussed later). 

The 77 and 77' mesons are often described in an SU(3) motivated quark basis, 
namely rys = {uu + dd-2ss)/V6, 771 = {uu + dd+ss)/V3. The mixing angle in 
this basis would be given by (j>— 54.7° in a lowest order chiral perturbation theory. 
In order to have fx ^ fn, one needs higher order terms in the chiral perturbation 
theory treatment and then the mixing scheme becomes more complicated Q| in 
this basis with more than one angle needed. 

In the SU(3) symmetric limit, m nn — m ss — m and x nn = x ns = x ss = x, 
so that only one mixing parameter is relevant and the mixing matrix simplifies 
considerably to a diagonal form with elements m 2 (octet) and m 2 + 3x (singlet). 
Previous lattice studies Q have used degenerate quarks, so have explored this 
case and have found that the mixing parameter x is of a magnitude which can 
explain qualitatively the observed splitting between the rj and rj mesons. 

Here we undertake a non-perturbative study in QCD from first principles 
which will be able to establish the values of the mixing parameters x, including 
the pattern of SU(3) breaking. This more comprehensive study would take into 
account the different masses of the light (u and d) quarks and the heavier s 
quark. Within the lattice approach, it is not at present feasible to evaluate 
using quarks as light as the nearly massless u and d quarks and also it is more 
tractable to use an even number of degenerate quarks in the vacuum. As we 
shall show, despite these restrictions, a thorough study of the mixing between 
7/ and 77' is possible. 

Our lattice study uses dynamical configurations with Nf = 2 flavours of sea 
quarks of type 1 and we consider the properties of pseudoscalar mesons made 
of either quark 1 or quark 2, where quark 2 corresponds to a heavier quark. 
Thus quark 2 is treated as partially quenched. Here we have in mind exploring 
a situation which will be relevant to treating strange quark propagation (quark 
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2) in a vacuum containing only lighter quarks (quark 1). We focus here on the 
results of lattice evaluations, for background to the methods used see rcf.Q. 
We address three topics where lattice input permits us to construct a firm 
foundation for the 77, rf mixing: 

• From comparing pseudoscalar meson masses with valence quarks of two 
different masses (namely meson masses mn, m\2 and TO22), we can esti- 
mate the mass m ss of the unmixed ss meson, given the observed m ns and 
ninn masses (ie K and tt respectively). 

• From measuring the mixing parameters X\\, X12 and X22 between initial 
and final flavour singlet states consisting of either quark 1 or 2 with dif- 
ferent masses as above, we can establish the pattern of SU(3) breaking in 
the mixing. 

• For Nj = 2 degenerate flavours of quark, we determine the pseudoscalar 
decay constants for the flavour singlet (Pq) and non-singlet (Pi) meson. 
This input allows us to discuss the relation between the observed 77 decay 
modes of 7r°, 77 and r( and the underlying quark content. 



2 Lattice results 

2.1 The ss pseudoscalar mass 

Chiral symmetry considerations lead to the expectation that the pseudoscalar 
meson composed of quarks of mass M q has mass squared m 2 which behaves 
linearly with M q at small quark mass. However, at large quark mass (c and b 
quarks for instance), one expects the meson mass to vary approximately linearly 
with the quark mass. Here we are not concerned with the region of very small 
quark mass where chiral logs are important jlj , so we summarise this behaviour 



For a pseudoscalar meson made of two different quarks of mass M n and M a , we 
shall assume its mass only depends on (M n + M s )/2 and not on (M s — M n )/2 
as found in lattice studies [Q] and in lowest order chiral perturbation theory. If 
eq. [3] were valid with just the linear term in the quark mass(ie c = 0), then 
one directly obtains the required mass of the pseudoscalar meson composed of 
s quarks, m 2 s = 2m 2 „ — m 2 m , that is 2if 2 — 7r 2 , leading to m ss = 0.687 GeV. 

This can be explored on a lattice by measuring the pseudoscalar meson mass 
for valence quarks in combinations 11, 22 and 12. Then, for small c/b, we have 



by 



m 2 = bM q + cM 2 q + 0{Ml) 



(3) 
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This has been studied in the quenched approximation giving evidence j5| for a 
positive coefficient c in eq. [j| In the quenched approximation, however, the chi- 
ral behaviour at small quark mass is anomalous since the theory is not unitary. 
A better way to study this issue on the lattice is then to use dynamical config- 
urations with sea quarks of type 1 and to consider the propagation of mesons 
made of either quark 1 or quark 2, where quark 2 corresponds to a heavier 
quark. 

We present results from UKQCD configurations || with Nf = 2 flavours of 
sea quark with SW-clover coefficient Csw — 1-76, lattice size 12 3 24, and with 
sea quarks having k = 0.1398, corresponding to sea quarks of mass around the 
strange quark mass {rap/ my = 0.67). Then we take the heavier valence quark 
(with k — 0.1380) as corresponding to approximately twice the strange mass 
(mp/mv = 0.81). 

The fits with two states to a 4 x 4 matrix of mesonic correlators for t range 
3 to 10 give results for the spectrum shown in Table |l|. 

Taking account of the correlation among the errors, we obtain the dimen- 
sionless ratio 

m 2 llC /(4b 2 ) = 0.011(3) (5) 

which indicates a statistically significant curvature from the c term. Setting the 
scale H using a -1 = 1.47 GeV then the value of c/b 2 in physical units can be 
obtained from ran = 698 MeV. 

Applying this value of c to the determination of the ra ss mass from the n 
and K masses, gives a relative shift upwards due to the curvature term (c) of 
1.1(3)%, corresponding to a value of m ss = 0.687 + 0.008 GeV. 

Note that this value also helps us to identify the meson mass ratio corre- 
sponding to strange quarks, namely rap/ray — iriss/ra^ = 0.682. 

This lattice study thus answers the question of the likely deviation in the 
pseudoscalar mass formula from the result given by the lowest order chiral ex- 
pression. 

2.2 Flavour-singlet mixing 

The mass splitting between flavour non-singlet and singlet mesons can be mea- 
sured using lattice evaluation of disconnected quark propagators. This is not an 
easy task: the contamination from excited states is difficult to remove and the 
statistical errors turn out to be relatively large. Initial studies have been in the 
quenched approximation ]2|, f^, 0]. Here, although there is no flavour splitting of 
the masses, the mass splitting matrix element x can be evaluated. It is, how- 
ever, preferable to be able to study the mass splitting directly and hence here 
we focus on results from full QCD simulations [|. 

The study of the mass spectrum of flavour singlet (Po) and non-singlet pseu- 
doscalar meson (Pi) using Nf = 2 flavours of sea quark 1 leads to singlet mass 
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0.1398 


0.1398 


0.1398 


0.477(5) 


0.171(6) 


0.56(4) 


0.196(12) 












0.56(5) 


0.176(18) 


0.1398 


0.1380 


0.1398 


0.563(5) 


0.182(4) 






0.1398 


0.1380 


0.1380 


0.640(4) 


0.190(4) 







Table 1: Pseudoscalar meson masses with Nf = 2, labelled 1 for non-singlet 
(isospin 1) and for singlet (isospin 0), and decay constants for the sea and 
valence quarks of hopping parameter shown (here 0.1398 corresponds to strange 
quarks and 0.138 to quarks twice as heavy as strange) from lattice studies with 
scale a -1 = 1.47 GeV. The singlet results are shown for fits with 1 state and 
3 < t < 7 (upper) and with 2 states and 2 < t < 7 (lower). 

mo = (^ii + 2IE11) 1 / 2 and non-singlet mass mj = rrin respectively which al- 
lows Xu to be extracted. We shall also be interested in the dependence of x on 
quark masses and on non-diagonal mixings. These can be studied with a little 
less rigour as we discuss later. 

We use the UKQCD lattices introduced in the previous section. The discon- 
nected diagrams were evaluated using a variance reduction method § which 
uses all the data available with no dilution from the stochastic method used. 
We use as many different operators for the pseudoscalar meson as possible to 
have the largest basis in which to extract the ground state - local and non-local 
(fuzzed) in space with both 75 and 7574 spin structure. Unfortunately, even 
with this basis of four operators, we are unable to determine the singlet mass 
precisely. For example for both valence and sea quarks having k — 0.1398, as 
shown in Table [j], we obtain am — 0.56(4) from a one state fit to t — 3 to 7 
with a 4 x 4 matrix of meson correlators. A two state fit to a wider t range 
(2-7) gives a similar mass value. The corresponding non-singlet pseudoscalar 
mass is also given in Table [lj so the determination of x from m\ + 2xu = m§ 
has relatively large errors (in = 0.10(4) GeV 2 using a^ 1 — 1.47 GeV). As an 
alternative, we also fit the ratio of the singlet to non-singlet correlators directly 
to a ground state mass difference. Using the t range 2-7 and local and fuzzed 
pseudoscalar operators, we obtain amg — ami = 0.12(7). This method gives a 
slightly larger mass value (indicating x\i = 0.13(8) GeV 2 ) but has even larger 
errors. These results indicate that much larger ensembles of gauge configura- 
tions will be needed to make more precise this approach of determining x from 
masses. 

If one studies correlations of meson operators made from valence quarks of 
type 2 in a sea of quarks of type 1, one will find the ground state pseudoscalar 
meson to be that composed of quarks of type 1 (assuming type 2 quarks are 
heavier than type 1). Because of this, we need to explore in more detail to study 
the SU(3) breaking of the mixing parameters. 
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X\2 


X22 


2 


.089(9) 


.073(6) 


.058(6) 


3 


.089(12) 


.072(10) 


.063(10) 



Table 2: Mass mixing matrix elements x in GeV 2 from quenched lattices for 
strange valence quarks (1) and quarks of twice the strange mass (2). 

To get a first look at this issue, we consider a quenched lattice and measure 
the ratio of the disconnected to connected diagrams for pseudoscalar meson 
propagation. We present results for (3 = 5.7, Csw — 1-57, 12 3 x 24 with 100 
configurations with k = 0.14077 and 0.13843. The non-singlet spectrum at these 
parameters was studied previously |^p|| giving my/mp values of 0.65 and 0.78 
which correspond approximately to strange quarks and quarks twice as heavy as 
strange. The scale was set as a" 1 = 1.2 GeV. The disconnected meson correlator 
was determined using a stochastic method with variance reduction j^] . 

Assuming dominance by ground state meson contributions, the ratio of dis- 
connected to connected diagrams at time separation t is 

£il = N fXij (t + l) 

dj 2(m ii m jj ) 1 / 2 [ ' 

with flavour non-singlet pseudoscalar mass ma for quarks of type i. The factor 
of t + 1 comes from the number of lattice sites at which the disconnected dia- 
gram can be split. To clarify the pattern of SU(3) breaking, we also study the 
non-diagonal case where we also measure the disconnected to connected ratio 
(here Cy is taken as (CuCjj) 1 ^ 2 ). In extracting X12 we can make an additional 
correction for the contribution from the propagation of mesons with different 
masses m\\ and 77122 although, in practice, this correction is very small. 

Using local meson operators, we obtain for x the values in Table |^. The 
largest t value has least contributions from excited state contamination and the 
consistency of the results versus t suggests that such contamination is small. As 
was found previously 0, x increases as the quark mass is decreased. Moreover, 
we can check to see if there is a factorisation of x as expected in some chiral 
perturbation theory descriptions JjJ, namely x\ 2 = ^11^22: and we find that X\ 2 
lies somewhat below the value given by this assumption. 

We now revert to discussing the more realistic (partially quenched) case: 
with heavier quarks of type 2 in a sea of two flavours of quarks of type 1 . The 
method described above for the quenched case can be applied here too. In 
principle this method is now only valid for small Nfxt for the propagation of 
quark 1. From this analysis of the measured D/C values, we get the x values 
shown in Table [| The values of x\\ are similar to those obtained above (with 
larger errors) directly from the rigorous method of using the mass differences. 
This suggests that the strong assumptions made in determining x directly from 
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.100(9) 
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.054(5) 
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.112(12) 
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.106(16) 


.077(15) 


.059(9) 


4F 


.093(13) 


.073(12) 


.052(11) 



Table 3: Mass mixing matrix elements x in GeV 2 from lattices with Np = 2 
flavours of sea quark (of type 1) with strange valence quarks (1) and valence 
quarks of twice the strange mass (2). The meson operators used are local except 
for the case labelled 4F which has non-local (fuzzed) creation and annihilation 
operators 

D/C are actually reasonable in practice. This, and the consistency of values 
from different t and different mesonic operators, gives us confidence to use the x 
values from quarks of type 2 (which are partially quenched anyway) as a guide 
to the quark mass dependence of x. The x values again show an increase with 
decreasing quark mass and also approximate factorisation. 

Setting the quark mass to strange (since mp/my = 0.682 in nature for s 
quarks) in both quenched and Nf = 2 evaluations leads to a consistent lattice 
estimate of x ss in the range 0.09 to 0.13 GeV 2 . This value is also consistent with 
that reported from a study of Nf = 2 by the CP-PACS collaboration jllj with 
mp/mv = 0.69 and a" 1 = 1.29 GeV giving values of x ss = 0.10 GeV 2 and 0.14 
GeV 2 (depending on using i m in = 2, 3 in fits, respectively). These lattice values 
are obtained at quite coarse lattice spacings and there may be some additional 
systematic error arising from the extrapolation to the continuum limit. We 
have, however, chosen to use a clover improved fermion action Q to minimise 
this extrapolation error. 

We are unable to determine the mixing strengths x for lighter quarks than 
strange. So we assume that the value of x continues to increase as the quark 
mass is decreased below strange in a similar way to the decrease we see from 
twice strange (type 2) to strange (type 1). 

Consider now the consequence of this determination of the mixing. We use 
input masses m nn — 0.137 GeV, m ss — 0.695 GeV (as discussed above) and 
aim to have x values in line with our results above, namely x ss ~ 0.12 GeV 2 , 
x\ s » XnnXss and we also expect, though with big errors from the extrapolation, 
Xnn/x ss ~ 2. The fi gure shows the x values needed to reproduce the known rj 
and rj masses for each mixing angle <f>. The lattice determination of x ss is shown 
by the dotted horizontal band. Keeping close to this band while satisfying the 
other lattice constraints is possible for the mixing illustrated by the vertical line. 
This has x nn — 0.292, x ns = 0.218, x ss — 0.13 GeV 2 which gives a description 
of the observed 77 and rf masses while being consistent with our QCD inspired 
evidence about the mixing strengths. This assignment corresponds to a mixing 
angle <f> in the r/ nn , i] ss basis of 44.5°. Note that this is almost maximal which 
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implies that the quark content (apart from the relative sign) of the rj and r\' 
meson is the same. The corresponding mixing angle in the r]g, fji basis (modulo 
comments above) is a value of 8 of —10.2°. 

2.3 Flavour-singlet decay constants 

The decays of 7r°, rj and if to 77 are expected to proceed via the quark triangle 
diagram. The quark model gives a decay proportional to Q 2 { for the contribution 
from a quark of charge Qi. Thus for the 7r° meson and the flavour-singlet nn 
and ss mesons, the quark charge contributions to the decay amplitudes would 
be in the ratio 1 : 5/3 : y/2/3. The experimental reduced decay amplitudes 
for 7T°, rj, and rj' are in the ratio 1.0 : 1.00(10) : 1.27(7). This information can 
be used to analyse the quark content of the pseudoscalar mesons subject to a 
quantitative understanding of the decay mechanisms. 

The conventional approach assumes that the decay constants for the decays 
of the three mesons are the same and then the relative decay amplitudes give 
information on the quark content. This suggests a mixing angle of 9 —20° is 
preferred % ||. 

We now address the issue of determining these decay constants directly from 
QCD using lattice methods. Our study uses 2 flavours of degenerate quark and 
we define the decay constants by 

<0|A"|Pi(<z)) - h<f (0|^o(5)) = f q* (7) 

For the isospin 1 state Pi (it - like), this is on a firm footing because of the 
anomaly hence f\ will be scale invariant. For the flavour singlet pseudoscalar 
meson the decay constant defined as above will not be scale invariant be- 
cause of gluonic contributions to the anomaly [Q . In this exploratory study we 
determine the decay constants with lattice regularisation and we shall compare 
the singlet and non-singlet values. 

These decay constants can be thought of as giving the quark wave function 
at the origin of the pseudoscalar meson. Since the mass splitting between singlet 
and non-singlet is not reproduced directly in quenched QCD, it is essential to 
use lattice studies that do include sea quark effects in this study of decay matrix 
elements. 

Results were obtained using fits to full (connected and disconnected) meson 
propagation with 4 different types of meson creation and destruction operator. 
These are local and fuzzed operators with either 75 or 7475 couplings, so giv- 
ing 4x4 matrix of pseudoscalar correlators. We used the Nf = 2 UKQCD 
configurations || referred to above. For the disconnected correlators, the vari- 
ance reduction technique § is essential to get a reasonable signal to noise ratio, 
particularly for the operators involving the 7475 factor. 

The lattice result for / with various valence quark masses with fixed sea 
quark mass (quark 1) as above is shown in Tabic [l| For the non-singlet results 
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using a -1 = 1.47 GeV and the tadpole-improved perturbative value of Z of 
0.81 (and of ca which is involved in mixing of the lattice pseudoscalar and axial 
currents but has a very small effect in practice) we get fn = 198(8) MeV. 
Since this corresponds to strange quarks, it is in reasonable agreement with 
experiment Q assuming a steady increase from /„„ = 131 MeV and f ns = 160 
MeV to f ss . We do see evidence for this increase in / with quark mass directly 
on the lattice going from quarks of type 1 (strange) to type 2 (twice strange) as 
shown in Table [|. 

The flavour singlet results are shown in Table [l]. They are determined by fits 
to the appropriate (connected plus disconnected) meson correlators which are a 
4x4 matrix at each t value. Despite this extensive data set, the determinations 
of / have relatively large statistical errors and the systematic error from chang- 
ing the type of fit is also comparable. For our case with Nf — 2 degenerate 
quarks, the comparison of the flavour singlet and non-singlet shows that the 
singlet decay constants appear to be somewhat larger, though the errors are too 
big to substantiate this. 

Combining the mass dependence we find in the flavour non-singlet sector 
with the near equality of singlet and non-singlet decay constants, we can de- 
duce properties of the physical case with three light quarks. Thus, in terms of 
the conventional treatment |l2|, we would expect > 1 and fr/'/fn > 1. 

One way to minimise the effects of mixing is to consider X = (a 2 + a 2 / )/ a 2 
where a refers to the reduced decay amplitude. Using the conventional formu- 
lae for the decay amplitudes would then give a value of X = 3r 2 (where r is 
a suitably weighted average of f v / and frf/fn which are both greater than 
1). Thus the conventional treatment gives X > 3 which is significantly larger 
than the experimental value |T^] of 2.64(24). Thus it appears unlikely that the 
conventional treatment (with the decay to 77 being given by the analogue of 
the formula for pions) is correct for any mixing angle. 

We conclude that there is no support for the conventional assumption that 
the singlet decays are given by a similar expression to the non-singlet. As has 
been pointed out by many authors jjj, this is plausible for at least two reasons: 
(i) the 77 and rf mesons are heavier and therefore less likely to dominate the 
axial current or, equivalently, higher order corrections to chiral perturbation 
theory will be more important (ii) the flavour-singlet axial anomaly has a gluonic 
component which will give additional contributions to any hadronic process. 

3 Conclusion 

From our careful non-perturbative study of mass formulae for flavour non-singlet 
pseudoscalar mesons made of different quarks, we deduce that the ss state lies at 
695 MeV. We then determine the pattern of mixing for the flavour singlet sector, 
obtaining x ss 0.12 GeV 2 , x nn /x ss w 2 and cc 2 s « x nn x ss . These conditions 
are indeed consistent and point to a mixing close to maximal (</> — 45 ± 2°) 
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in the nn, ss basis (this corresponds to a conventional (77s, 771) mixing 9 of 
— 10 ± 2°). We are able to explore the decay constants for singlet pseudoscalar 
mesons for the first time. Our results show similar decay constants for singlet 
and non-singlet states of the same mass but with quite large errors. 

We have not addressed here the issue of the origin of these mixing parameters 
x. Lattice studies || have the capability to relate them to topological charge 
density fluctuations or to other vacuum properties. 

Our lattice studies have been hampered by two constraints. One is that the 
disconnected quark diagrams needed for a study of singlet mesons are intrinsi- 
cally noisy. Much larger data sets (tens of thousands of gauge configurations) 
will be needed to increase precision. Another constraint is that we are unable 
to work with sea quarks substantially lighter than strange. We have also not 
attempted a continuum limit extrapolation of our lattice results. Although we 
are using a lattice formalism that should improve this extrapolation, it would 
be safer to test it directly. The lattice non-perturbative results do, however, 
show clearly the structure of the mixing in the singlet pseudoscalar mesons. 
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Figure 1: The mass mixing parameters x in GeV 2 versus 77, 77' mixing angle <f> in 
the i]nn, Vss basis. The horizontal dotted lines give the allowed range from the 
lattice determination of x ss . The vertical line illustrates our preferred solution. 
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